function plot_bellows_drift_correction_toolbox(patient, fig_num)
    
     %| Plot original vs drifted bellows voltage, using the patient.drift
     %| as;
     %|     patient.bellows_volt_drifted = patient.bellows_volt + patient.drift*patient.scan_times_zerod; 
    
    %------------------------------------------------------------------------
    %|      Dependancies;  
    %|                  
    %|                 
    %|                 
    %------------------------------------------------------------------------
    %   This file is part of the
    %   5D-Novel4DCT Toolbox  ("Novel4DCT-Toolbox")
    %   DH Thomas, Ph.D
    %   University of California, Los Angeles
    %   Contact: mailto:dhthomas@mednet.ucla.edu
    %------------------------------------------------------------------------
    % $Author: DHThomas $	$Date: 2014/04/01 10:23:59 $	$Revision: 0.2 $
    
    if nargin < 2
        figurecount = length(findobj('Type','figure'));
        fig_num = figurecount+1;
      
    end
    
    figure(fig_num);
        set(gcf,'units','normalized','position',patient.default_figure_position);

    hold on;
    p1 = plot(patient.scan_times_zerod(:),patient.bellows_volt(:),'.');
    p2 = plot(patient.scan_times_zerod(:),patient.bellows_volt_drifted(:),'r+');
    hold off
    legend([p1(1,1);p2(1,1)],'Original',sprintf('Drift Corrected = %.2f [mV/s]',patient.drift*1E3),'Location','NW')
    xlabel('Scan Time [s]')
    ylabel('Bellows Volt [V]')
